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Few interacting fermions in one-dimensional harmonic trap 
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We study spin-1/2 fermions, interacting via a two-body contact potential, in a one-dimensional 
harmonic trap. Applying exact diagonalization, we investigate the behavior at finite interaction 
strength, and discuss the role of a ground state degeneracy which occurs for sufficiently strong 
repulsive interaction. Even low temperature or a completely depolarizing channel may then dra- 
matically influence the system's behavior. We calculate level occupation numbers as signatures of 
thermalization, and we discuss the mechanisms to break the degeneracy. 



Introduction. Advances in the control and manip- 
ulation of ultracold quantum gases have opened up a 
new avenue to the study of interacting few particle sys- 
tems [1]. In particular, existing trapping techniques aUow 
for exploring the physics in low dimensions where the 
quantum-statistical distinction between fermionic and 
bosonic particles experiences severe modifications. A 
striking property of one-dimensional systems is that a 
mapping between a strongly repulsive bosonic system to 
a non- interacting fermionic system exists |2|-|j|, giving 
rise to a strongly correlated phase known as the Tonks- 
Girardeau (TG) gas. In a gas of ultracold *^Rb atoms, 
this "fermionization" has for the first time been observed 
experimentally [5]. Also the strongly attractive coun- 
terpart of the TG gas, a strongly correlated metastable 
phase known as the super TG gas [6|, has been observed 
recently with ^■^'^Cs atoms |7|. 

Similar physics is also exhibited by spin-1/2 fermions in 
a one-dimensional trap: A fermionic analog of the super 
TG phase has been discussed in Ref. [^]. In the strongly 
repulsive limit, the spin-1/2 fermions may form a ground 
state which is identical to the one of non-interacting 
fermions without spin |9i] . High-precision control of such 
systems has been proven feasible in a recent experiment 
which allows for preparing the system in a state with a 
well-defined, small number of particles |10| . In particu- 
lar, it has become possible to study the ground-state and 
dynamics of a two-fermion system [11], for which the ex- 
act theoretical solution is known in the full interaction 



parameter range |12l.ll3|. 

For a description of systems with three or more 
fermions, different analytical and numerical methods 
have been applied [14,-1 18.], suggesting such systems as 
a tool for studying ferromagnetism, and providing some 
insight in the fermionized nature of the strongly repulsive 
system. A crucial difference between the fermionization 
of a spin- 1/2 Fermi system and the TG fermionization 
of the spinless Bose gas is the fact that the fermionic 
Girardeau state is only one candidate amongst many de- 
generate states in the limit of infinitely strong repulsion. 
How the system behaves in this limit thus remains un- 



clear: Will there be a thermal state, as temperature be- 
comes large compared to the small degeneracy splitting 
for large but finite interactions? Is there a completely 
depolarizing channel which according to the Jaynes prin- 
ciple would demand an equal superpositioir of all possible 
ground states? 

With available techniques, it will be possible to di- 
rectly address these questions experimentally [10]. Here 
we give a theoretical description, based on an exact di- 
agonalization study, discussing possible signatures of the 
degeneracy in an experiment. In particular, we show that 
different superpositions of the ground states lead to sig- 
nificantly different occupation numbers of the harmonic 
oscillator levels. These quantities, accessible to experi- 
ment via an opening of the trap, are shown to distin- 
guish betweeir thermal states at different temperatures, 
and the maximum-eirtropy state. 

A closely related question concerns the mechanisms 
which could provide the mixing of the degenerate states: 
As shown by our numerics, starting from a unique ground 
state at finite repulsion, both an adiabatic increase or a 
quench of interactions will bring the system into only one 
state of each degenerate manifold, which is protected by 
permutation-group symmetry against mixing with other 
states from the quasi-degenerate manifold. For an un- 
derstanding of the system behavior, it thus turns out 
to be relevant to consider mechanisms beyond the ideal 
model which are able to break the degeneracy: First we 
show that anharmonicities in the trap have no infiuence, 
as they are independent from the spin symmetry of the 
state. We then explicitly show that a Zeeman effect due 
to an inhomogeneous magnetic field is able to mix dif- 
ferent spin symmetry sectors, giving rise to a degener- 
acy lifting. A similar effect can also be caused by spin- 
dependent p-wave interactions and spin-orbit coupling. 

System. Our system consists of two-species fermions 
confined in a one-dimensional trap. The Hamiltonian 
has the form 
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where gm is an effective interaction strength between 
two fermions of different spins, and V = ^mui'^x'^ is a 
harmonic trapping potential. Here we assume that both 
species have equal masses m, and are confined by the 
same trap with frequency w. We will thus, in the follow- 
ing, refer to the two species as a single species with an 
internal (pseudo) spin- 1/2 degree of freedom. For further 
convenience, energies will be expressed in units of huj, 
and lengths in units of the harmonic oscillator length 

Xosc = y^h/mUJ. 

Let us note that the interaction term in the Hamilto- 
nian ([T]) is non-zero only for states with a spatial wave 
function which is symmetric under particle exchange. For 
fermions with the same spin the wave function is always 
antisymmetric, and the interaction term will not con- 
tribute. For two fermions with opposite spins symmetric 
and antisymmetric wave functions are possible, and cor- 
respond to states with zero and finite interaction energy. 
This observation will play a crucial role in the following 
part of the paper. 

The single-particle problem yields a basis (t>n{x) for 
our studies of the many-body problem. These functions 
are simply the harmonic oscillator eigenfunctions corre- 
sponding to a energy e„ = n -I- 1/2. We introduce opera- 
tors a„CT and a\„^ which annihilate and create a fermion 
with spin a in the level n. The Hamiltonian then reads 



if = ^ e„c 
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where for convenience we introduce g — guj/huj, and 
fmnip are the matrix elements of the contact interaction 
in the harmonic oscillator basis. The properties of the 
Hamiltonian ([T]) are studied with the exact diagonaliza- 
tion method. This can be done by fixing the total num- 
ber of particles N, and the numbers N^ (-^4.) of t ii) 
fermions, and truncating the single-particle basis at some 
sufficiently large level rimax (we take rimax = 20). 

Before turning to our numerical results, let us con- 
sider two limiting cases which can be solved analyti- 
cally. The first case is the trivial one of g = 0, that 
is, a non-interacting system. It is clear that the ground 
state in this case is obtained by simply filling the lowest 
rmn{N-^,Ni) levels with two fermions of opposite spin. 
The remaining fermions will singly occupy the subse- 
quent levels. Without losing generality we set the energy 
of this state E-p — 0. 

The second limiting case is the Girardeau limit of 
infinitely strong repulsive interaction between the two 
species, g — s> cxd. Then, a Fermi-Fermi mapping, in anal- 
ogy to the Bose-Fermi mapping introduced by E. Lieb 
and W. Liniger [J|, allows one to treat the repul- 
sive two-species fermions like non-interacting one-species 
fermions. A spatial wave function for the ground state 
is then obtained as a Slater determinant of the N lowest 
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FIG. 1: We plot the energy E — Ef as a function of the 
interaction strength g for diflferent combinations of spin-up 
and spin-down particles. The energy offset _Ep is the Fermi 
energy of the non-interacting system. We find ground state 
degeneracies in the limit of strong interactions. We plot first 
20 eigenenergies of the Hamiltonian. 



levels. It can be rewritten as [9|: 
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This spatial wave function is fully antisymmetric, and 
thus corresponds to a fully symmetric spin configuration. 
It is an eigenfunction of both the single-particle and the 
interaction part of the Hamiltonian, and thus provides 
an exact eigenfunction for any choice of 5. In particu- 
lar, as the wave function vanishes whenever two particles 
are at the same position, it describes a state with an en- 
ergy which is independent from g, and which becomes 
the ground state energy for g -^ 00. 

It is possible to symmetrize Eq. ([3]) with respect to 
pairs of particles of opposite spin just by including a fac- 
tor sgn{xi — Xk)- However, as our numerical results sug- 
gest, wave functions obtained in that way are eigenfunc- 
tions of the Hamiltonian only for g — ;> 00. More insight 
is provided by an exact solution of the two-particle prob- 
lem ^2,]. This is done by separating the center-of-mass 
part of the Hamiltonian. The relative of motion of the 
two particles is then described by 



1 d- 
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^rcl = -7r3-2+7rr^ + <?<5M 



2dr2 
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where r — ri — r2 is the relative position of two par- 
ticles. The relative motion part of the wave function 
^ is found to be the first excited state of the Hamil- 
tonian Q with energy E = 1 for any g. Its center- 
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FIG. 2: (a) The probabilities Pi, P^ for finding t and I par- 
ticles in different orbitals (n — 0, f , 2, 3) as a function of the 
interaction strength g. We consider the effect of temperature 
using a Boltzman distribution: the solid black and dashed 
blue lines denote the temperatures kT/huj = and 0.3 respec- 
tively, (b) The probabilities Pj_ of finding 4- particle in given 
level as a function of temperature for g = 12. All probabilities 
saturate in the limit of completely depolarizing channel. In 
both plots N = 5 and A''4, = 1. 



of-mass motion is in the ground state. In the limit 
g — >■ cxD, the state ^ becomes degenerate with the 
ground state which smoothly evolves to the symmetric 
state ^o(''7^) oc |r|e~'' /^e~^ /^, as we adiabatically in- 
crease g. Note that this wave function, despite describing 
a state of zero interaction energy, is not an eigenstate of 
Eq. (ID) for any finite g. 

We thus have seen that in the limit of infinitely strong 
interactions, the two-particle problem has two degener- 
ate ground states with opposite symmetry of the spatial 
wave function. One solution is obtained from the other by 
multiplication by sgn(xi — 12). This operation turns the 
spatially antisymmetric wave function ^ into a spatially 
symmetric wave function, and thus has to be accompa- 
nied with a corresponding change in the symmetry of the 
spin wave function. 

The same mechanism can be applied for larger systems, 
N > 2. Then, for every pair of particles with opposite 
spin, it is possible to change the symmetry of the spatial 
wavefunction in the state ^ , and thereby construct new 
solutions in the Girardeau limit. This has been done in 
Ref. [14;], and leads to a degenerate ground state mani- 
fold, where the number of degenerate ground states D is 
given by the number of distinct spin configurations. It is 
counted by the distinct possibilities of dividing N parti- 
cles into two groups with A^-|- and N^ members, that is, 
D = jyT^f I ■ Note that the degeneracy of higher man- 
ifolds (related to the excitations of the center of mass 
of the system) changes since one also has to take into 
account excitations in the relative motion. 



Role of degeneracies. These degeneracies, although 
known before |l4| . have, to our knowledge, not been 
studied. As we show below, they might, however, play 
a crucial role in understanding the few-body physics of 
strongly interacting fermions. 

Our numerics focuses on the region between the two 
limiting cases, where, for N > 2, exact solutions are not 
known. In that region, the system makes use of both 
the possibility of doubly occupying the lowest levels to 
reduce potential energy, and occupying higher levels in 
order to reduce interactions. Energy spectra as a function 
of interaction strength g are plotted in Fig. [TJ 

We find different energy manifolds which become de- 
generate in the limit <; — >■ 00. Each manifold corresponds 
to different center-of-mass wave functions. The number 
of degenerate states in the lowest manifold is given by D. 
Contrarily, at any finite g, this degeneracy is lifted. The 
highest energy state of the lowest manifold is described 
by the fully antisymmetric wave function of Eq. ^ . As 
explained above, it is an exact solution with zero inter- 
action energy for any g. Its energy function is therefore 
simply a horizontal line. 

We next consider the population of the different single- 
particle levels. In the quasi-degenerate regime, the van- 
ishing small energy gap does not protect the ground state 
against mixing with other states from the manifold: If 
the gap and the system's temperature are of similar or- 
der, the system might best be described by a thermal 
state. 

As we show in the upper line of Fig. [5] for N^ — 4 
and iV^ = 1, a finite temperature then strongly affects 
the occupation probabilities of the 4, particle. Taking 
into account the whole manifold of five quasi-degenerate 
states, we apply a Boltzman average at different tem- 
peratures. For g > 10, the probability of finding the 4, 
particle in the level n = is clearly reduced even by 
a small temperature k^T <^ hu. Similar effects can be 
found in other correlation functions, such as, for instance, 
the joint probability of finding an ti pair iu one of the 
harmonic-oscillator levels. 

To understand this, we note, that at zero temperature, 
as shown by the thick solid line, only the true ground 
state is occupied. For this state, the probability for the 
4, particle to be in n = ranges between 1 at 5 = down 
to 0.8 for g ^ 1. In contrast, the Girardeau state has 
an equal population of the five lowest level, such that the 
probability of finding the 4- particle in n = is given by 
0.2. This shows that the quasi-degeneracy of the ground 
state enhances the population of the higher energy levels 
in the limit 17 ^- 00. 

Interestingly, as shown by the middle line of Fig. [2l 
the occupation probabilities for the t particles are close 
to unity, independently from temperature and interaction 
strength g, as a consequence of the Pauli principle. 

As shown in Fig. |31 the energy gap within the quasi- 
degenerate manifold decays exponentially fast as a func- 
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FIG. 3: (color online) Energy gap of the quasi-degenerate 
manifold as a function of the interaction for different number 
of particles and different imbalances. The black solid line is 
the energy gap between the ground state and the first excited 
state. The red dashed line is the width of the ground mani- 
fold. Both quantities decay exponentially with the interaction 
strength g. 



tion of g. The thermal regime can thus be expected to 
quickly approach the scenario in which the temperature is 
infinitely large compared to the vanishingly small many- 
body gap, but still small on the single-particle energy 
scale huj. In this limit we can assume the existence of 
a completely depolarizing channel favoring a state with 
maximum entropy. As shown in the lower line of Fig. [51 
for an intermediate interaction strength, the completely 
depolarizing channel leads to occupation numbers which 
are clearly distinguishable from those at finite tempera- 
ture. 

External symmetry breaking. So far we have assumed 
thermalization mechanisms that are able to bring the 
system into a superposition of different states from the 
quasi-degenerate manifold. On the basis of our ideal 
model, however, it turns out that neither an adiabatic 
increase nor a sudden quench of the interaction param- 
eter lead to occupation of more than one state in the 
degenerate manifold. This is the state with the largest 
total spin, and, in accordance with the Lieb-Mattis the- 
orem [19|, the state with the smallest energy. A mixing 
with other states is prohibited by the permutation-group 
symmetry, that is, the conservation of total spin. 

Not surprisingly, as shown in Fig. IH^a), a possible an- 
harmonicity in the trap only shifts the quasi-degenerate 
energy manifolds towards larger energies without lifting 
the degeneracy. From the nature of the symmetry, only 
mechanisms which simultaneously act on spin and spatial 
degrees of freedom may break the degeneracy: This could 
be a spin-orbit coupling, spin-dependent interactions like 
for instance p-wave interactions, or the existence of a 
spatially dependent magnetic field. The latter option is 
easily implemented as a Zeeman term Hz — S^^Xiaf 
in the Hamiltonian, where we, without loss of generality, 
assume that the two internal states have opposite mag- 
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FIG. 4: (color online) (a) Influence of a symmetric anhar- 
monicity of the trap to the spectrum of the Hamiltonian: 
Anharmonicity shifts eigenenergies of the Hamiltonian (black 
lines) relative to the eigenenergies of the Hamiltonian with 
harmonic trap (grey lines). It does not lift the degeneracies 
at large g, and provides no mixing between different states in 
each manifold, (b) Energies as a function of g in the presence 
of a Zeeman term: No degeracies occur in the strongly inter- 
acting limit, (c) For A'^ = 4 and A'^^ = 1, the ground state 
at 3 = 12 is time-evolved after switching on a Zeeman term, 
(5 = 0.05. The probability for finding the state in one of the 
four possible spin configurations is plotted. The thick black, 
dashed blue, dotted red, and thin black lines correspond to 
the states ordered by energy (increasing). This illustrates the 
Zeeman breaking of the spin symmetry. 



netic moments along z-direction, and (5 is a magnetic field 
gradient. Furthermore, it can also be implemented in a 
controlled way in the experiment. We will thus discuss 
the effect of a Zeeman term in more detail. 

First, we consider two particles with opposite spin. It 
is easily seen that the symmetric wave function |S) ex 
|r|e~'' /^e~^ 1"^ (which has to be multiplied by an anti- 
symmetric spin wave function \\\^ — \\Xl) has non-zero 
transition matrix elements (S|iJz|A) — 2-\/2/7r(5 with 
the antisymmetric wave function |A) ex re"*" '^e~^ " 
(which has to be multiplied by a symmetric spin wave 
function |ti) + lit))- This gives rise to a degeneracy 
splitting A which is linear in the magnetic field gradient, 
A = 4-y/2/7r(5, and to mixed-symmetry states in the limit 
of large g. 

The same effect is found for systems with three or four 
particles, where we have taken into account the Zeeman 
term in our exact diagonalization study. The situation is 
plotted in Fig. Sfb), clearly showing the lifted degener- 
acy in the large g limit. In particular, we find that for 



sufficiently small Zeeman energy, 6{x) <i; fiuj, only states 
of the same energy manifold are mixed. Furthermore, 
since the Zeeman term can be rewritten as a sum of op- 
erators acting only on pairs of particles, transition matrix 
elements of states which differ by more than one unit of 
total spin are zero. Accordingly, the matrix representa- 
tion of the Zeeman term has a tridiagonal structure. 

The mixed symmetry of the eigenstates in the pres- 
ence of a Zeeman splitting is illustrated by Fig. IDJc): 
For TV = 4 and N^ = 1 we consider a system which is 
prepared in the (maximum total spin) ground state of the 
Hamiltonian for g = 12 and S — 0. Then we switch on 
the external magnetic field gradient 6 = 0.05, and propa- 
gate the state for some time t. We then measure the spin 
symmetry of the state (by projecting back into the spin- 
conserving basis given by the eigenstates of 5 = 0). We 
then plot the probabilities of finding the system in one 
of the three different spin sectors as a functions of time. 
The time scale of the dynamics shown in Fig. \M,c) can 
be controlled by the strength of the field gradient. This 
could allow for studying the crossover from quantum time 
evolution to thermalization. 

Conclusions. In our study of one-dimensionally 
trapped spin-1/2 fcrmions we have focussed on the 
strongly repulsive regime, in which a ground-state de- 
generacy is exhibited. This quasi-degenerate manifold 
allows one to study thermalization in a small quantum- 
mechanical system. We have calculated correlation func- 
tions and occupation numbers of the harmonic oscillator 
levels as experimentally accessible signatures for distin- 
guishing between pure states and thermal states at finite 
or even infinite temperature. Since each eigenstate is 
protected from mixing with other states by its symme- 
try with respect to the permutation group, mechanisms 
for thermalization must in general be operators which 
simultaneously act on spin and spatial degrees of free- 
dom. We have shown that the presence of an additional 
Zeeman term lifts the degeneracy, and may lead to a 
time-dependent superposition of different states from the 
quasi-degenerate manifold. 

Note. When finishing this work, we became aware of a 
related work by S. E. Gharashi and D. Blume exploring 
the degenerate regime of strongly repulsive ID fermions 
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